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I.  Introduction 


The  modeling  of  surface  acoustic  wave  transducers  and  reflectors  has 
been  discussed  by  several  authors.  Each  model  is  restricted  to  transducers 
or  reflectors  and  most  are  valid  over  a  limited  frequency  range.  This  work 
provides  energy  conservation  based  corrections  for  a  mixed  units  scattering 
matrix  [1,2]  recently  derived  from  an  infinite  array  analysis.  This  thesis 
is  the  first  publication  to  present  results  involving  transducers  with 
floating  electrodes  and  is  a  general  analysis  which  models  any  periodic 
device  over  all  frequencies  of  interest. 

This  thesis  presents  a  universal  analysis  of  surface  acoustic  wave 
transducers  and  finite  length  reflectors.  The  final  result  of  this  one 
dimensional  analysis  is  a  3  x  3  mixed  units  scattering  matrix  description 
of  a  reflector  or  transducer  whose  parameters  are  functions  of  frequency. 

This  analysis: 

*  is  applicable  to  all  periodic  structures,  where  metallization 
ratio  and/or  electrode  spacing  may  vary  slowly  in  the  direction 
of  propagation; 

*  models  interdigital  transducers,  regardless  of  connection 
circuitry,  including  arbitrary  voltage  electrodes,  floating 
electrodes,  and  electrically  loaded  reflector  arrays; 

*  is  valid  for  all  frequencies  of  interest; 

*  is  based  on  analytical  expressions  involving  material 
constants,  electrode  geometry  and  individual  electrode  circuitry; 

*  provides  the  input  admittance,  electrical  to  acoustic  transfer 
function  and  acoustic  to  acoustic  transfer  function  (transmission 
and  reflection); 

*  is  in  a  form  that  directly  accepts  a  recently  derived  set  of 
mechanical  scattering  coefficients  [27,28]. 
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A  number  of  useful  techniques  have  been  developed  to  analyze  SAW  reflectors 
and  transducers.  A  brief  summary  of  these  techniques  and  their  features 
adds  insight  to  the  capabilities  of  this  mixed  units  scattering  matrix 
approach . 

Realizing  that  the  location  of  interdigital  electrodes  corresponds  to 
the  time  delayed  signal  generated  by  an  incident  acoustic  impulse,  several 
authors  have  described  an  impulse  response  model  [ 3 » 8 » 29  3 .  This  model  is 
particularly  applicable  to  weak  coupling  piezoelectric  materials  where  tap 
to  tap  interactions  are  neglected.  The  electro-acoustic  transfer  function 
of  each  tap  is  assumed  to  be  flat  over  the  bandwidth  of  interest  and  the 
tap  weights  are  the  Fourier  transform  of  a  desired  frequency  response, 
making  first  order  modeling  a  simple  process.  Hartmann  extended  this  model 
to  include  the  calculation  of  transducer  impedances  and  the  effect  of 
electrical  loading  on  filter  response  (second  order  regenerated 
reflections)  via  an  electrical  port  equivalent  circuit  [4,6].  He  points 
out  that  another  frequency  dependent  acoustic  source  could  be  Engan's 
Fourier  series  solution  [7]  for  the  electrostatic  field  produced  by  the 
electrodes  of  an  infinite  array,  but  limited  his  analysis  to  the 
fundamental  response.  Arrays  which  use  overlap  weighting  with  a  fixed 
voltage  are  modeled  using  a  channelized  approach  similar  to  one  described 
by  Tancrell  and  Holland  [8]. 

In  their  initial  strong  coupling,  transmission  line  analysis  Smith  et 
al.  [4]  presented  two  infinite  array  models  ("in-line"  and  "crossed  field") 
derived  from  the  Mason  bulk  wave  transducer  equivalent  circuit  [5].  The 
selection  between  the  two  approaches  is  based  either  upon  the  predominance 
of  an  electric  field  normal  or  parallel  to  the  surface.  Smith  and  Pedler 
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[9]  extended  this  work  by  introducing  1)  quarter  wavelength  transmission 
line  sections  of  different  impedances  (values  are  empirically  derived) 
representing  the  electrodes  and  gaps  to  model  acoustic  reflections,  2) 
Engan's  source  field  into  their  crossed  field  Mason  bulk  wave  model  to 
analyze  harmonics,  and  3)  missing  electrodes  and  end  effects  cases  to 
handle  finite  length  arrays.  Their  source  function  is  the  derivative  of 
the  normal  electric  field  component  with  respect  to  the  direction  of 
propagation  along  the  surface.  This  new  source  field  extends  the  analysis 
accurately  to  odd  harmonic  operation.  The  field  equations  were  solved  for 
a  number  of  electrode  connection  sequences,  taking  into  account  next 
nearest  neighbors.  The  results  are  presented  in  tabular  form  and  are  handy 
for  design  of  periodic  polarity  sequence  transducers.  The  authors  suggest 
that  the  case  of  floating  electrodes  with  induced  charges  is  handled  by 
numerically  resolving  the  electrostatic  boundary  value  problem. 

Although  Smith  and  Pedler's  tables  are  accurate  near  odd  harmonics, 
Szabo  realized  the  need  for  inter-harmonic  frequency  analysis  and  applied  a 
Fourier  transform  of  the  displacement  field  as  the  source  function  in  the 
spectral  weighting  model  [10].  A  set  of  electro-acoustic  transfer 
functions  are  calculated  numerically  using  field  theory  for  a  number  of 
local  electrode  environments.  These  frequency  domain  curves  replace  the 
flat  passband  response  of  the  delta  functions  used  in  the  impulse  response 
models  and  serve  as  the  mathematical  link  between  such  models  and 
experimental  device  behavior.  The  curves  are  used  in  a  design  prescription 
to  calculate  the  tap  weights  for  a  given  frequency  response.  Closed  form 
expressions  of  these  transfer  functions  for  structures  other  than  arrays  of 
single  electrodes  are  difficult  to  obtain  analytically. 
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The  advent  of  floating  electrodes  and  capacity  weighted  transducers 
brought  about  the  need  for  a  new  model  which  does  not  require  numerically 
recalculating  the  field  solution  for  each  possible  electrode  voltage 
configuration.  Datta  et  al.  [11]  point  out  that  in  general,  for  the 
impulse  response  model,  the  separation  of  the  transducer  response  into  an 
element  factor  and  an  array  factor  is  not  possible,  since  the  charge 
distribution  on  each  electrode  is  different  because  of  neighboring 
electrode  effects.  They  calculated  an  element  factor,  referred  to  as  the 
basic  charge  distribution,  defined  as  the  charge  induced  in  a  grounded 
array  with  one  volt  applied  to  the  center  tap.  The  basic  charge 
distribution  does  not  change  from  tap  to  tap  and  extends  the  impulse  model 
to  all  harmonics,  taking  next  nearest  neighbors  into  account.  This 
approach,  based  on  closed- form  field  theory  solutions,  accurately  models 
the  electro-acoustic  response  of  arbitrary  electrode  voltages  and  apodized 
transducers.  The  weak  coupling  approximation,  on  which  this  work  is  based, 
does  not  allow  the  prediction  of  acoustic  reflections  and  frequency  shift 
due  to  the  local  electrode  shorting  of  the  acoustic  potential  wave  on 
highly  piezoelectric  substrates. 

Aoki  and  Ingebrigsten  [12]  recognized  the  importance  of  reflections 
and  frequency  shift  in  an  infinite  array.  They  derived  the  equivalent 
circuit  parameters  of  the  "mixed  circuit"  model  (a  modification  of  Smith's 
"in-line"  model)  from  the  closed  form  field  theory  dispersion  relations 
near  resonance  to  obtain  a  non-eopirically  based  transmission  line  model. 
Although  neighbor  coupling  effects  when  arbitrary  voltages  are  applied  (as 
in  floating  electrodes  and  capacity  weighted  transducers)  are  neglected, 
this  dispersion  relation  based  (strong  coupling)  model  correctly  predicts 
the  transducer  acoustic  reflection  coefficient  magnitude  as  a  function  of 


frequency.  This  analysis  agrees  well  with  experiment  near  the  odd 
stopbands.  The  cascaded  "mixed  circuits"  do  not  produce  a  stopband  at  the 
even  harmonics  as  field  theory  suggests,  therefore  Aoki  and  Hattori  [13] 
added  a  stored  energy  "shunt-susceptance  transmission- line"  model  for  use 
at  the  even  harmonics  and  showed  good  agreement  with  experiment  for  special 
transducer  configurations. 


A  traveling  acoustic  wave  creates  a  set  of  phased  voltages  on  an  array 
of  electrodes  if  acoustic  reflections  are  assumed  negligible.  The  phase 
difference  between  electrodes  is  constant  for  a  given  frequency.  The 
scattering  matrix  of  a  single  electrode  in  this  infinite  periodic,  phased 
voltage,  environment  has  been  calculated  from  field  theory  by  Datta  and 
Hunsinger  [1],  The  scattering  matrix  elements  are  analytical  expressions 
which  are  functions  of  frequency,  metallization  ratio,  and  substrate 
parameters.  The  analysis  works  well  on  reflectors  and  multistrip  couplers 
(passive),  but  does  not  address  the  case  of  transducers  in  which  an 
arbitrary  voltage  is  applied  to  each  electrode,  nor  to  fixed  tap  weight 
transducers  at  non-synchronous  frequencies. 


These  single  electrode  piezoelectric  scattering  coefficients  for 
infinite  periodic  arrays  [ 1 ]  have  been  modified  for  modeling  finite  length 
periodic  transducer  and  reflector  structures.  A  single  electrode  is 
represented  as  a  3  x  3  mixed  units  scattering  matrix.  These  matrices  are 
cascaded  (via  acoustic  ports)  to  form  an  N+2  by  3  matrix  that  describes  the 
device.  This  matrix  is  then  compacted  into  a  3  x  3  matrix  which  completely 
defines  the  three  port  transducer  or  two  port  finite  length  reflector 
without  the  use  of  special  case  look-up  tables  or  curves.  The  analysis  is 
based  on  piezoelectric  and  mechanical  parameters,  accounting  for  both 
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scattering  mechanisms  in  slowly  varying,  periodic  structures  with  arbitrary 
voltages. 

In  section  II.  the  modifications  to  the  single  element  scattering 
parameters  along  with  the  conditions  required  for  cascading  are  described, 
while  the  cascading  procedure  is  described  in  section  III.  The 
interpretation  of  the  composite  scattering  matrix  is  discussed  in  section 
IV.  along  with  a  detailed  comparison  of  the  analysis  with  experimental 
results  for  reflectors  and  single  (solid)  and  double  (split)  electrode 
transducers  with  varying  loads  at  any  frequency.  This  section  includes  a 
comparison  of  experiment  with  theory  when  the  added  mechanical  scattering 
is  taken  into  account.  Transducers  with  floating  electrodes  are  discussed 
in  section  V.,  with  a  conclusion  in  section  VI. 


II.  Single  Electrode  Mixed  Units  Scattering  Parameters 

This  section  describes  how  the  scattering  parameters  [a']  for  a  phased 
array  transducer  (Fig.  la)  are  applied  to*  determine  the  mixed  units 
scattering  matrix  [a]  for  a  single  electrode  unit  cell  in  a  grounded  array 
(Fig.  1b).  The  mixed  units  scattering  matrix  characterizes  the  three  port 
unit  cell.  The  two  acoustic  ports  are  characterized  by  positive  and 
negative  traveling  acoustic  waves,  while  the  electrical  port  is  defined  in 
terms  of  the  most  conveniently  measured  parameters,  voltage  and  current. 

The  mixed  units  scattering  parameters  for  the  single  electrode,  unit  cell 
are  of  the  form 


where  A^  is  the  positive  propagating  wave  (open  surface  electrical 

t*  H 

potential)  as  it  enters  the  i  electrode,  Bi+1  is  the  negative  propagating 
wave  and  1^  is  the  current  induced.  The  unit  cell  has  the  properties  of 
being  symmetrical  and  lossless.  The  parameters  m^  s  n^,  m^  =  m^  ,  and 
®12  =  ^23  are  dimensionless  while  m^  =  m^2  and  m^  have  the  units  of  mhos. 
Also,  m^  is  related  to  m^  by  a  constant  (derived  in  Appendix  I.B.). 

The  single  electrode  scattering  matrix  is  evaluated  from  the 
previously  derived  phased  array  results  by  noting  that  the  single  tap  and 
phased  array  transducers  are  identical  when  the  impressed  voltage  (Vg) 
equals  zero.  This  allows  the  first  two  columns  of  the  [m]  matrix  to  be  set 
equal  to  their  phased  array  matrix  counterparts  (equation  (23)  in  [1]). 


When  a  voltage  is  applied  to  one  electrode  in  a  grounded  array  a 
charge  distribution  is  induced  under  all  adjaoent  electrodes.  Acoustic 
waves  are  generated  by  each  of  the  charged  electrodes  and  the  total  wave 
generated  is  the  sum  of  all  the  contributions  (A')  as  shown  in  Fig.  1c. 

The  electro-acoustic  transfer  function  (a^  a  m^)  of  the  ith  electrode  in 
a  grounded  array  is  defined  in  teras  of  the  total  acoustic  wave  generated 
due  to  the  voltage  applied.  The  model  assumes  that  the  entire  wave  is 
generated  in  the  ith  electrode.  Each  contribution  of  A'  is  phase  shifted 
to  the  center  of  the  1th  electrode  and  suaaed  to  give  the  total 
electro-acoustic  transfer  function.  The  equations  for  a^  and  mg^  nr*  the 
charge  distribution  (o  (s,A  ),  equation  (7)  in  [2])  aultipied  by  the 
factor  (  A v/v)(p/^}+6Q)  where  A  v/v  is  half  the  piezoelectric  coupling 
coefficient,  p  is  the  spacing  between  electrodes  and  €p  +€q  is  the 
effective  substrate  dielectric  constant. 

The  current  entering  the  electrode  when  no  incident  acoustic  waves  are 
present  is  represented  by  a  complex  number  (m^,  derived  in  section  II.B.). 
The  real  part  arises  due  to  electrical  to  acoustic  power  conversion  which 
is  proportional  to  the  coupling  coefficient.  The  imaginary  part  is  the 
electrostatic  capacitive  admittance  based  on  the  non-piezoelectric 
properties  of  the  material.  For  simple  periodic  voltage  transducers 
(ie.  single  or  double  electrode)  this  capacitance  is  the  capacitance  per 
pair  (for  example,  see  equation  (14)  in  [6]),  while  for  devices  with 
arbitrary  voltage  electrodes,  an  inter-electrode  capacitive  admittance 
matrix  (derived  from  equation  (8)  and  (9)  in  [2])  is  added  during  the 
cascading  process  (section  V.). 
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The  week  coupling  scattering  parameters  Cm']  (referenced  to  the 
electrode  center)  are  tabulated  here  for  convenience. 


Acoustic  dependent 


o' 
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a  m' 
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Voltage  dependent 


a' 
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(2d) 


where 

Pn  is  the  nth  Legendre  polynomial 
Av/v  is  the  half  the  piezoelectric  coupling  coefficient 
N  *  integer (2f/fQ) 

fQ  a  open  surface  velocity/ electrode  spacing  »  v/2p 
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a  »  2f/fo  _  N 

A  *  tt  Tl  (  T]  being  the  metallisation  ratio) 

a  is  the  normalized  radian  frequency 

g  is  the  substrate  dielectric  constant 
P 

W  is  the  acoustic  beam  width 

The  a,,,  values  are  found  from  equations  (9)  and  (24)  in  [1].  The  derivation 
°f  is  found  in  section  II.  B. 

This  [m*]  matrix  uses  the  weak  coupling  approximation  which  assumes 
that  the  incident  (or  generated)  wave  amplitude  is  constant  while 
traversing  the  electrode  region.  This  approximation  is  quite  good,  however 
a  close  inspection  of  the  [m1 ]  matrix  will  show  that  energy  conservation  is 
not  precisely  maintained  (jm'^l2  ♦  !m* 21  i 2  *1).  Any  errors  in  the 
original  single  electrode  matrix  are  accumulated  and  magnified  by  the 
cascading  process.  Therefore,  the  initial  terms  [m*]  are  corrected  (to 
[m])  by  accounting  for  wave  growth  under  the  electrode.  A  first  order 
effect  (2b)  is  the  phase  shift  (5  =  jA)  of  m'2r  As  will  be  shown,  the 
second  order  correction  requires  a  decrease  in  the  magnitude  of  m'gj.  The 
problem  separates  into  the  shorted  case  (acoustic  terms  only)  and  the 
externally  electrical  loaded  electrode  case. 

1 

l 
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A.  Shorted  Case 

The  shorted  case  involves  only  the  acoustic  terms  of  the  mixed  unit 
scattering  matrix  [m]  (essentially  an  acoustic  2-port),  therefore  only  the 
upper  left  corner  is  addressed 


The  conditions  of  losslessness  for  a  2-port  require  that  the  net  input 
power  to  a  lossless-passive  junction  be  zero  for  any  combination  of 
excitation  possible.  That  is,  there  is  no  energy  absorbed  in  the  Junction. 
Therefore,  all  of  the  terms  of  the  2x2  scattering  hermitian  matrix  are 
set  to  zero  (equation  (2.22)  in  [23]).  Hence,  there  are  four  lossless 
conditions: 


(4a) 


(4b) 


(4c) 
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whdre  zoi  311(1  Z02  are  the  transmission  line  Impedances  at  their  respective 
ports. 

These  equations  with  Zq^  equal  to  Zg2  for  a  symmetrical  two  port  (m^ 

8  m22  311(1  ml2  *  m21^  reduce  to  two  equivalent  and  more  revealing  forms  to 
be  satisfied  (equation  (2.25)  in  [23]): 

1  -  jm1 t i 2  -  !m21 1  2  =  0  (5a) 

and 

®12  *  9n  +  (2n-1)  n/2.  (5b) 

The  modification  of  the  acoustic  terms  involves  calculating  corrected 
values  of  1)  the  magnitude  of  m21,  2)  the  angle  of  m21  and  3)  the  angle  of 
m^.  The  transmission  vector  of  a  unit  cell  without  a  scattering  element 
(no  electrode)  is  the  unit  vector  t'  in  Fig.  2.  When  an  electrode  is 
introduced  into  the  unit  cell,  part  of  the  transmitted  wave  is  continuously 
scattered  as  the  wave  travels  under  the  electrode.  The  electrode  is 
conceptually  divided  into  many  point  sources,  each  being  a  scatterer  of 
acoustic  waves  (small  arrows  in  Fig.  2a).  The  reflections  from  each 
section  (reverse  scattering)  experience  phase  shifts  (with  refernce  to  the 
electrode  center)  and  therefore  are  summed  as  phasors  (Fig.  2a)  to  produce 
The  weak  coupling  approximation  assumes  that  the  transmitted  wave 
traverses  each  section  of  the  electrode  unscattered  and  accounts  for  the 
total  forward  scattering  with  the  addition  of  an  imaginary  scattered  wave 
(JA)  generated  by  the  electrode  (Fig.  2a).  This  term  is  negative  due  to 
the  localized  electric  field  shorting  which  slows  the  wave  as  it  propagates 
under  the  electrode.  When  6  =  JA  is  added  to  the  unit  transmission  vector, 
the  energy  balance  (5a)  is  not  satisfied  because  the  magnitude  of  m'21  (= 


i  1 


Figure  2.  Scattering  coefficient  vector  addition,  acoustic  terms, 
shorted  case,  (a)  Weak  coupling  approximation,  (b)  Energy  conservation 
correction. 


.5 


!l  ♦  j 6  1 )  is  greater  than  one.  Therefore,  6  is  rotated  (see  Fig.  2b)  to 
form  one  side  of  a  triangle  composed  of  'Si,  the  unit  vector  and  the 
magnitude  of  m^1  (as  yet  unknown). 

The  first  correction  is  the  calculation  of  the  transmitted  wave  vector 
magnitude.  This  second  order  correction  (decrease  in  magnitude  of 
®'2-|)  is  very  small.  The  magnitude  of  acoustic  reflection  (m'^)  is  also 
affected  by  this  second  order  correction  (ie.  change  in  all  scattering 
element  magnitudes),  but  since  lm'^1  <<  Im'^l  the  absolute  error  in 
la' 11  I  is  near  zero  (as  is  the  error  for  all  ImVjl,  ij  +  12  or  21).  The 
magnitude  of  m21  is  then  calculated  using  (5a)  and  the  magnitude  of  m'.^. 

The  second  step  involves  calculating  the  phase  angle  of  the 
transmitted  wave  vector  angle  (921)  using  the  unit  transmission  vector, 

!  6  !  and  I m2 -j  !  -  The  law  of  cosines  is  applied  to  the  triangle  formed  by 
these  three  vectors  (Fig.  2b) 


giving  the  angle  of  referenced  to  the  electrode  center.  The 
transmission  phase  delay  across  the  unit  cell  is  implemented  by  adding  a 
tt  f/fQ  phase  shift,  thus  giving  the  corrected  angle  of  m2r  Equation  (5b) 
for  a  symmetrical  two  port  results  in  m2i  and  m22i.  m^.  Therefore  for 

the  third  correction  (5b),  the  angle  of  m^  is  adjusted  to  be  tt/2  radians 
from  921<  The  two  port  acoustic  scattering  matrix  now  satisfies  the 
lossless  conditions. 


.  k  V  I* 
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B.  Externally  Electrical  Loaded  Electrode  Case 


At  this  point  there  are  5  corrections  remaining:  the  angles  of  a^, 
a^,  m^., ,  a^  and  the  .complex  value  of  m^-  Talcing  the  solution  of  m^ 
first,  from  energy  considerations  it  is  shown  in  appendix  I. A.  that  the 
electrical  port  radiation  conductance  is 

Ga  s  Refa^)  =  --5(21^0*31  +  ffl23n*32^  ^ 


If*  Z.  is  defined  as  the  characteristic  electro-acoustic  impedance: 


z0  =  -i-_LL 

2  P. 


where  A  i3  the  open  surface  electric  potential  and  P  is  the  total 
*  ® 

acoustic  power  flow,  then  using  (appendix  I.B.) 

-,3-  - 


we  find  that 


Ga  ,  2  mi 3m*3l 


for  a  symmetrical  lossless  3-port  and  is  independent  of  any  scattering 
coefficient  angle.  The  imaginary  part  of  electrical  port  input  admittance 
arises  from  both  electrostatic  (capacitive)  and  acoustic  generation.  The 
acoustic  contribution  (8  )  is  the  Hilbert  transform  of  the  real  electrical 

cl 

port  admittance  [30].  For  a  single  electrode  the  value  of  B  is  near  zero 

cl 

because  the  Hilbert  transform  of  a  slowly  varying  G  is  very  small. 

cl 

Therefore,  the  imaginary  part  of  m33  is  only  the  admittance  due  to  the 
capacitance  of  the  unit  cell  electrode.  For  a  single  electrode  transducer 
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the  capacitance  of  an  electrode  pair  (Cp)  is  equal  to  the  capacitance  of 
one  electrode  in  an  alternately  connected  array.  The  total  single 
electrode  admittance  is  then 

O33  =  2  m1 3°*31  ♦  J<“Cp  (11) 

Z0 

The  solution  of  the  4  unknown  electro-acoustic  angles  is  now 
addressed .  The  problem  is  simplified  by  the  property  of  a  symmetrical 
3-port  in  which 

“3.,  =  m^2  (12a) 

and 

®13  =  “23  (12b) 


These  two  equations  along  with  (9)  reduce  the  problem  to  one  unknown  angle 
(for  instance,  9^)  fro®  which  the  three  other  angles  are  found. 


The  solution  of  an  equation  containing  9^  (r  arg(m13))  involves 
finding  the  2-port  matrix  equivalent  to  the  3-port  matrix  when  the  third 
port  is  electrically  terminated  in  a  load  (*L).  Mathematically,  the  load 
is  absorbed  into  the  third  port  electrical  admittance.  Setting  the 
electrical  port  external  current  equal  to  the  load  current  (-VY^)  and 
solving  for  the  normalized  voltage  at  that  port  induced  by  an  incoming 


(13) 


The  last  equation 
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-VYl  =  b31a1  ♦  033V 


(14a) 


after  re-arrangement  and  the  use  of  (8),  produces  the  normalized  voltage 
V  _  -m^  _  -2  m1 

A, 


11 


'l  m33+YL  Z0(m33+YL) 


(4b) 


Inserting  this  relation  into  the  first  row  as  the  generating  voltage  gives 


the  reflection  coefficient  (m^  g^) 


1  =  mii  *  mi3v  =  mn 


m.,V  _  m.,  -  2m 
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:  m 


Z0  ^m33+YL^ 


1  lei 


(15a) 


Similarly,  for  the  transmitted  wave 


*2  =  “2 


•  2m 


II 


=  a. 


Z0  ^m33+YL^ 


'2  lei 


(15b) 


Therefore  the  electrically  loaded  acoustic  scattering  matrix  (supression  of 
the  third  port)  is  represented  as 

r a,i 

(16) 


—  — 

r-  1 

— 

B1 

_ 

m1  lel 

m12el 

A1 

A2 

m21el 

m22el 

* 

MOT  __ 

L  J 

L  J 

The  same  symmetrical  lossless  2-port  equations  (5a, b)  apply  to  this 
matrix.  Equation  (5a) 

“"W2  *  '"2lel'2  -  '  *  0  <17a) 


is  now 


„  2 

2 

„  2 

«„  -  2  ’’3 

m  2  1,, 

m21  “  13 

Z0  (.33.^) 

Z0  (a33+rL) 

-1=0  (17b) 
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which  simplifies  to 

1 2  m  ^  I 

ReCm^  +.  m21)  cos  013  +  Im(m1  ^  +  m21)  sin  0^  I  13  [  -  0 

!Z0  (m33+YL)| 

(17c) 

giving 

I 

-1 


ei3  =  2 


0,-,  +  0 


33  *  (m1 1+m21 ) 


-  cos 


12  ®13  /(Z0(n>33+YL))| 


'm1 1  +  °21 


(17d) 


The  resulting  angle  and  the  use  of  (9)  and  (12)  give  the  values  of  the  four 
unknown  angles.  As  a  check,  the  angle  of  °-]iel  is  compared  with  the  angle 
of  m2iel  311,1  f°un<*  t0  satisfy  (5b). 

In  summary,  the  corrected  mixed  units  scattering  coefficients 
referenced  to  the  edges  of  the  unit  cell  (Fig.  1b)  are: 


”13  = 
n31  = 


m22  s  I m *  1 , 1  e^912^ff/2)  e^^o 
m12  =  (  1  -  !m* , ,  ! 2)  1/2  eJ012  e_Jllf/fc 
m23  =  !m'  i  eJ013  e’^^o 


m32  =  l»’31!  eJ(013  -n)  e^/2fo 


(18a) 

(18b) 

(18c) 

(18d) 


°33  =  Ga  +  J"  Cp 


(18e) 


These  single  electrode  scattering  matrix  coefficients  conserve  energy 
for  acoustic  reflection,  generation  and  detection.  In  the  calculation  of  a 
finite  length  device  matrix,  many  single  electrode  matrices  are  cascaded 
and  the  final  device  matrix  satisfies  the  three  port  energy  conservation 
equations  (equation  (2.165)  in  [26]). 
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III.  Composite  Transducer  Matrix 

A.  Cascading 

In  order  to  analyze  an  entire  transducer,  the  effects  of  individual 
fingers  and  their  interactions  must  be  combined,  while  taking  into  account 
the  finger  connection  sequence.  This  is  done  by  cascading  3-port  mixed 
units,  single  electrode  [m]  scattering  matrices  ( 1 8 )  and  then  compacting 
the  final  matrix  into  a  3  x  3  matrix  [M]  describing  the  total  array 
(Fig.  3). 


"B1 " 

M11  M12  M13 

■*r 

an+i 

3 

M21  M22  M23 

0N+1 

w 

_ 1 

_M31  M3  2  M33_ 

_VT_ 

The  scattering  matrix  cascading  process  starts  with  a  single  element 


scattering  matrix  [m*-]  (Fig.  4,  equation  (20)).  This  matrix  is  combined 


with  another  adjacent 

Cm2] 

matrix 

which  may 

be 

different 

than 

the  firs 

_  i 

1 

1  _ 

_  2 

2 

2  -| 

B1 

®11 

®12 

“13 

"  A1~ 

"V 

®11 

”12 

“l3 

"A2~ 

A2 

- 

1 

°21 

®22 

1 

®23 

B2 

A3 

- 

2 

®21 

2 

®22 

2 

®23 

B3 

_1 

Jj31 

®32 

®33_ 

_h_ 

_X2_ 

2 

J?31 

2 

®32 

2 

®33L 

_V2_ 

(20) 


The  acoustic  interface  between  the  two  unit  cells  contains  two  unknown 
interacting  acoustic  waves  A2  and  B2>  The  two  sets  of  scattering  equations 
(20)  also  have  two  equations  relating  A^  and  B2.  These  2  equations  and  2 
unknowns  are  solved  for  Ag  and  B^  in  terms  of  A^  B^,  V1  and  V2-  These  are 


j,  •'  >.-v  ■* 


Figure  3.  Generalized  weighted  transducer  driven  by  balanced  voltage 


source . 
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substituted  into  their  respective  locations  and  a  new  4x4  scatter  matrix 
is  produced  which  relates  the  electrode  pair  input  quantities  to  the  output 
quantities. 


m’n 

m'l2 

m*  1 3 

m'  14 

CM 

B 

m'22 

“’23 

“*24 

m’31 

m*32 

m'33 

a*  34 

m*  4 1 

a'42 

m*  43 

“*44 

A1 

B3 

V1 

V2 

(21) 


This  new  matrix  describing  the  electrode  pair  now  becomes  the  left 
matrix  and  another  3x3  matrix  representing  the  third  electrode  is 
combined  with  it.  Each  matrix  has  an  equation  relating  the  interacting 
acoustic  potentials  (now  and  B^)  in  terms  of  ,  B^,  V^,  V2,  and  V^. 

The  2  equations  with  2  unknowns  are  solved  and  a  new  5x5  scattering 
matrix  for  the  electrode  triplet  is  created.  This  process  continues  until 
the  final  N+2  by  N+2  array  scattering  matrix  is  found  for  an  N  electrode 
device  (Fig.  3)- 


"Bi  - 

m11  m12  n13  “l4  *  *  * 

A1 

AN+1 

m21  °22  m23  n24  ’ 

BN+1 

h 

“31  ®32  |  “33  m34 

V1 

CM 

M 

“41  “42  |  “43  “44  *  *  ‘ 

V2 

• 

•  •  j  •  •  •  •  • 

• 

• 

•  •  j  •  •  •  •  • 

• 

• 

•  •  |  «  •  ••• 

• 

•  •  |  •  •  •  •  • 

_VN_ 

I 
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As  stated  previously,  for  devices  with  periodic  electrode  voltages,  the 
capacitive  admittance  is  included  in  the  single  electrode  scattering 
matrix.  For  devices  with  arbitrary  voltages,  a  sub  matrix  of  capacitive 
admittances  (Dr  from  equations  (8),  (9)  in  [2])  is  added  to  the  admittance 
sub  matrix  due  to  acoustic  generation  (dashed  line  in  (22)).  That  is,  for 
i>2  and  j>2 

mij  total  8  Bij  acoustic  + 

This  matrix  (22)  completely  describes  the  interactions  between  the  incoming 
acoustic  traveling  wave  potentials,  the  voltages  impressed  on  the 
electrodes  and  the  outgoing  acoustic  traveling  wave  potentials  and  the 
electric  port  induced  currents.  Appendix  II.  contains  a  more  detailed 
description  of  the  cascading  process. 
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B.  Compacting 

The  array  scattering  matrix  can  describe  any  periodic  surface  wave 
device  taking  into  account  the  connection  sequence.  In  an  interdigital 
transducer  (IDT)  the  electrode  connection  sequence  (and  therefore  each  tap 
voltage)  is  known  apriori.  Each  of  the  electrode  voltages  is  related  to 
the  terminal  voltage  (V^)  by  a  voltage  tap  weight  a^  (Fig.  3)»  where 

ais(W  (24) 

Therefore  in  (22),  m1  ^>2  may  be  replaced  by 
1 3  8  m13a1  +  °l4a2  +  ®15a3  + 

(25) 

m,23  s  m23a1  +  “24a2  +  ®25a3  * 


and  the  n>2  columns  are  compacted  into  one. 
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• 
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Therefore,  if  the  tap  voltages  are  known,  the  largest  matrix  is  N+2  by  3 
(as  shown  in  appendix  II.). 

This  matrix  is  compacted  by  row  as  well.  The  terminal  current  (IT  in 
Fig.  3)  is  the  sum  of  the  currents  in  the  electrodes  connected  to  one  bus 
bar.  This  current  compacting  is  approached  in  much  the  same  way  as  the 
voltage  compacting.  Let  b^  be  the  current  tap  weight,  that  is  the  portion 
of  the  electrode  current  generated  by  the  incident  acoustic  wave  that  flows 
into  the  positive  bus  bar  (ie.  high  side,  1^  in  Fig.  3).  In  most  IDT 
structures  (but  not  capacitive  weighted)  b^  will  be  either  0  or  1 .  In  a 
single  electrode  IDT  every  other  b^^  is  1  because  every  other  electrode  is 
connected  to  the  high  side.  Therefore,  i>2  is  rePlaced  by 

®"31  a  +  °4l^2  +  ®5l^3  +  *** 

m"32  =  m32b1  +  m42b2  +  m52b3  +  . . .  (27) 

®"33  s  m'33b1  +  m'43b2  +  m'33b3  *  ... 

and  the  m>2  rows  are  compacted  into  one. 


— 1 
w 

®1 1 

®12 

®'  13 

AN+1 

- 

®21 

®22 

®'23 

_Jt  _ 

■"31 

®"32 

“"3i 

N+1 


(28) 


The  result  is  the  desired  3x3  mixed  units  scattering  matrix  (19)  which 
completely  describes  the  transducer  as  a  three  port. 


.  L ...  ‘lanlMl  -  ’■ 


j. 
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IV .  Scattering  Matrix  Analysis  of  Transducers  and  Reflectors 

The  knowledge  of  the  total  scattering  matrix  for  a  given  surface  wave 
device  tap  configuration  is  a  very  powerful  analysis  tool.  A  number  of 
examples  will  be  examined  to  show  its  usefulness. 


A.  Transducer  Acoustic  Reflection  Coefficient 

The  acoustic  reflection  coefficient  of  an  electrically  loaded 
transducer  is  derived  from  the  array  scatter  matrix  (19),  (28).  The 
shorted  transducer  (V^  -  o)  reflection  coefficient  due  to  MEL  reflections 
is  M11#  The  reflection  coefficient  of  an  electrode  (or  transducer)  loaded 
by  an  admittance  is  composed  of  two  parts,  the  MEL  reflection  and  a 
regenerated  wave.  This  regenerated  wave  is  produced  by  the  induced  voltage 
at  the  electrical  port  arising  from  an  incident  wave .  To  calculate  this 
voltage  for  the  case  of  a  transducer,  BN+1  is  set  equal  to  zero,  IT  is  the 
load  current  (-V^Y^)  and  VT  is  solved  for  as  in  the  single  electrode  case 
(13),  O1*)  and  (15).  The  M  matrix  equations  are  written  as 

B1  =  M  A  *  MV 

'd  i  (29) 

"vTYL  =  M31A1  +  ^M33+*L^VT 

Solving  the  second  equation  for  the  terminal  voltage  gives 

VT  =  ~M31  A1 

M33+YL  (3°> 


Inserting  this  into  the  first  equation  results  in 


B1  ■  M11  *i  -  m13m31  *1 
"33**1. 


(31a) 


* 

Reflections  from  shorted  electrodes  (Mechanical-Electrical  Loading). 
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and  the  electrically  loaded  reflection  coefficient  is 


rel  =  M11  "  m13M31 


M33+YL 


(31b) 


As  an  example,  the  complex  reflectance  for  a  20-pair,  single  electrode 
(Fig.  5,  inset),  0.5  duty  factor  (metallization  ratio)  uniform  weighted 
transducer  on  Y-Z  lithium  niobate  [14]  is  plotted  for  two  loads  in  Fig.  5. 
These  calculations  show  a  calculated  frequency  dependent  reflection 
coefficient  whose  shape  has  good  agreement  with  experiment.  The  authors 

[15]  attributed  a  similar  difference  in  magnitude  between  theory  and 
experiment  to  the  measurement  method.  This  calculated  response  also  agrees 
with  that  calculated  by  Aoki  and  Ingebrigsten  (Fig.  4  in  [12]). 

B.  Transducer  Radiation  Conductance 

The  transducer  radiation  conductance  Ga  (the  real  part  of  M^) 
includes  the  effects  of  regeneration,  stored  energy  and  other  tap 
interactions.  The  input  conductance  of  a  78-pair,  0.5  duty  factor, 
unweighted  linear  transducer  on  Y-Z  lithium  niobate  has  been  calculated  and 
is  compared  in  Fig.  6  with  the  experimental  results  for  a  30  MHz  device 

[16] .  The  theoretical  calculation  agrees  very  well  with  experiment. 
Previous  theories  did  not  predict  the  velocity  change  under  the  electrodes. 
Although  Jones  et  al.  had  calculated  the  magnitude  accurately,  they  were 
not  able  to  predict  the  1.7$  frequency  shift  down  from  f  without  any 
post-calculation  fitting. 


(ftW)  (D0)BO| 
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Calculated  radiation  conductance  for  a  15-pair,  0.5  duty  factor, 
double  electrode  transducer  is  plotted  in  Fig.  7  and  shows  good  agreement 
with  experiment  [20].  The  center  frequency  of  this  device  is  assumed  to  be 
49.44  MHz. 


C.  Reflector  Analysis 


The  scattering  matrix  approach  lends  itself  readily  to  the  analysis  of 
reflectors.  The  acoustic  reflection  coefficient  of  a  shorted  reflector  is 
the  same  as  that  of  a  shorted  transducer.  The  reflection  coefficient  of  an 
open  isolated  electrode  reflector  is  unlike  that  of  an  open  circuited  (no 
load)  transducer,  although  the  reflection  mechanisms  (MEL  and  regenerated 
wave)  are  the  same.  In  the  case  of  isolated  electrodes,  each  electrode 
operates  as  an  independent  (no  electrical  load)  reflector  whose  Induced 
voltage  is  a  function  of  the  incident  acoustic  potential.  It  is  conve¬ 
nient  here  to  introduce  a  neutral  reference  potential.  For  an  IDT  this 
neutral  reference  is  located  half  way  between  electrodes  of  opposite  polarity, 
so  that  one  half  the  terminal  voltage  exists  between  each  terminal  and  neu¬ 
tral.  Admittances  are  now  defined  from  each  electrode  to  the  neutral  refer¬ 
ence.  At  fundamental  the  capacitance  of  each  electrode  to  the  neutral 
reference  is  twice  the  electrode  to  electrode  capacitance  used  previously. 


“33  iso 


*  2J*Cp 


(32) 


The  single  electrode  three  port  is  reduced  to  a  two  port  by  supressing  the 
electrical  port.  The  electrically  loaded  equations  (14-16)  (Y^  =  0)  are 
evaluated  with  the  isolated  admittance  (32)  and  replace  the  shorted 
acoustic  terms  (3)  in  the  single  electrode  scattering  matrix  [m] 
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This  two  port  matrix  completely  defines  the  open  circuit,  three  port,  unit 
cell  in  an  isolated  array. 


As  an  example,  consider  the  case  of  one  electrode  in  an  infinite 
array.  At  center  frequency  the  single  electrode  acoustic  reflection 
coefficient  magnitude  for  the  open  and  short  circuit  cases  are  equal  (at 
0.5  duty  factor)  and  on  Y-Z  lithium  niobate  the  value  of  0.01647  agrees 
well  with  experimental  results  [16-19].  For  the  case  of  double  electrodes 
(ie.  0.5  duty  factor  electrodes  connected  in  pairs),  the  calculated 
coupling  coefficient  at  fundamental  is  reduced  by  a  factor  of  0.736  from 
the  single  electrode  case.  This  predicts  a  piezoelectric  transducer 
coupling  factor  K  =  (0. 046X0. 736) 1 =  0.184  compared  to  the  experimental 
value  of  0.193  [12]. 

To  find  the  reflection  and  transmission  coefficients  of  isolated 
electrode  arrays,  the  single  electrode  scattering  matrix  (33)  is  cascaded 
as  in  section  III. A.  (the  non-acoustic  terms  are  set  equal  to  zero).  The 
calculated  array  acoustic  transmission  coefficients  (shorted  and  open 
cases)  are  compared  with  experiment  [21]  (bars)  for  three  duty  factors  in 
Fig.  8.  The  device  is  a  100  electrode  array  on  Y-Z  lithium  niobate.  The 
experimental  values  of  single  electrode  normalized  impedance  (z  in  [21]) 
are  used  in  the  author's  transmission  line  equation  (with  N  s  100  and  duty 
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Figure  8.  Transmission  coefficient  versus  duty  factor  for  a  shorted 
and  open  (isolated)  100  electrode  reflector  on  Y-Z  lithium  niobate. 
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factor  s  t  (author's  notation)  =  0.5)  to  calculate  the  array  reflection 
and  transmission  coefficients: 


R  =  2'  -  1 
Z'  ♦  1 


IT!  =  (1  -  |R!2)1/2 


where  z'  =  z2*1  for  the  short  circuit  and  z'  =  z”2N  for  the  open  circuit 
case.  Good  agreement  with  experiment  is  found  in  Fig.  9  where  the 
transmission  coefficient  (T  =  lO-log^M^ ))  of  this  100  electrode  array  is 
plotted  versus  frequency. 

The  device  transmission  coefficient  (Mg^)  is  a  complex  number  whose  phase 
describes  the  dispersion  of  finite  length  reflector  and  transducer  arrays. 

The  fractional  acoustic  transmission  coefficient  phase  shift  (^21  ^21 ) 
defined  as 

4*21  *21 


Ifnf/f. 


where  Nrrf/fo  is  the  open  surface  phase  delay.  This  fractional  phase  shift 
is  equal  to  the  fractional  shift  in  velocity  iv/v  and  is  compared  with 
experiment  (duty  factor  =  1)  =0.5)  [22]  in  Fig.  10  for  an  80  electrode 
device  on  Y-Z  lithium  niobate  at  three  duty  factors.  This  calculation  is 
for  three  possible  electrode  connection  sequences:  all  connected  (shorted) 
in  Fig.  10a,  every  other  electrode  connected  together  (open  circuit  IDT)  in 
Fig.  10b,  and  all  electrodes  unconnected  (isolated  array)  in  Fig.  10c.  The 
highly  dispersive  nature  of  the  curves  near  the  fundamental  and  second 
harmonic  is  evidence  of  a  large  number  of  reflections  in  the  arrays.  This 
region  is  denoted  as  a  stopband.  It  is  interesting  to  note  that  the 


Figure  10.  Normalized  relative  acoustic  transmission  phase  shift 
versus  frequency  (80  electrodes)  for  three  duty  factors,  experiment  (x  x  x) 
at  0.5  duty  factor,  (a)  Shorted  array,  (b)  Open  circuit  transducer,  (c) 
Isolated  array. 


shorted  and  open  transducers  (Fig.  10a  and  b,  respectively)  have  identical 
dispersion  characteristics  out  of  this  band.  The  value  of  &v/v  at  the 
fundamental  stopband  corresponds  to  the  1.7$  fractional  shift  in  frequency 
evident  in  all  the  devices  discussed  (for  example,  see  Fig.  6).  All  three 
connection  sequences  display  the  presence  of  a  second  harmonic  stopband 
which  is  also  shifted  downward  in  frequency. 

D.  Mechanically  Loaded  Transducers 

A  recently  derived  method  of  calculating  the  mechanical  scattering 
coefficients  [27,28]  allows  this  scattering  matrix  cascading  process  to 
include  the  effects  of  elastic  loading,  mass  loading,  and  stored  energy  in 
thick  electrodes,  grooves  and  buried  electrodes.  The  mechanical 
coefficients  (as  in  the  piezoelectric  case  (2a,b))  are  reflection 
coefficient  and  transmission  phase  shift  for  the  unit  cell  and  are  composed 
of  contributions  due  to  linear  (impedance  mismatch)  and  quadratic  (stored 
energy)  effects 

"11  oech  *  ♦  r2<h/»2) 

2  (36) 

J*21  mech  s  J<¥hA>  *  Vh/X)  > 

and  at  the  present  are  valid  for  frequencies  about  the  harmonics.  The 
derivation  of  these  constants  (r1(  r g,  Ag)  is  found  in  reference  28 
(specifically  equations  24  and  26)  and  first  order  examples  of  the 
calculations  are  found  in  reference  27.  The  mechanical  coefficients  (36) 
are  added  to  the  initial  piezoelectric  scattering  coefficients  (from 


(2a, b),  where  JA  21 

piezo 

s  jA  from  m'21 

)  to  give  the  composite 

terms 

“’ll 

total 

8  “'ll  piezo 
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which  are  used  in  the  energy  balance  (5,6),  thereby  effecting  a  complete 
scattering  analysis. 

An  example  of  the  mechanical  scattering  contribution  is  shown  in 

Fig.  11  where  the  calculated  input  radiation  conductance  is  plotted  versus 

frequency.  This  is  compared  with  experiment  [24]  for  a  10-pair,  single 

electrode  interdigital  transducer  (duty  factor  s  0.5).  The  major 

reflection  mechanism  of  the  transducer  in  Fig.  11a  with  1500  ?  thick 

aluminum  electrodes  is  piezoelectric  loading  (localized  shorting  of  the 

traveling  potenial  wave  under  the  electrode) .  The  effect  of  combined 

o 

piezoelectric  and  mechanical  loading  is  shown  in  Fig.  11b  where  750  A  of 
gold  has  been  added  to  the  transducer  of  Fig.  11a.  Although  the  calculated 
magnitudes  are  high,  the  frequency  shifts  of  the  upper  zeroes  are 
accurately  predicted  by  this  analysis.  Many  authors  have  attributed 
lowered  experimental  magnitudes  to  scattering  into  bulk  waves  and  viscous 
dampening  of  the  surface  wave. 


(a)  1500  A  thick  aluminum  electrodes. 


V.  Transducers  with  Floating  Electrodes 


Transducers  whose  electrodes  are  individually  connected  to  external 

loads  and/or  each  other  are  analyzed  in  a  straight-forward  manner  using 

this  generalized  scattering  matrix  approach.  The  discussion  of  section 

III.  dealt  with  electrodes  whose  voltages  were  known  apriori,  whereas  this 

section  deals  with  transducers  which  contain  contain  floating  (unconnected) 

electrodes  in  combination  with  connected  electrodes.  Reflector  arrays  of 

floating  (isolated)  electrodes  have  been  treated  in  section  IV. D.  where  the 

current  through  each  electrically  unloaded  electrode  was  a  known  value 

(zero).  This  section  addresses  itself  to  a  similar  problem  (each  floating 

electrode  current  equals  zero) ,  but  uses  circuit  analysis  techniques  to 

determine  the  unknown  voltages.  These  voltages  are  then  used  to  derive  the 

3x3  matrix  [M]  by  solving  for  the  individual  scattering  parameters  in  a 

manner  similar  to  the  compacting  sequence  of  section  III.  Although  the 

floating  electrodes  discussed  are  electrically  isolated,  this  is  in  fact  a 

specific  case  of  electrodes  which  are  externally  electrically  loaded.  In 

another  case,  the  electrodes  could  be  connected  in  groups  to  a  load  (load  cu 

rent  *  -V_YT)  or  to  individual  loads  Y.  ,  where  the  load  currents  equal  -V  Y 
T  L  i-i  1 

As  an  example,  consider  a  solid  electrode  interdigital  transducer  with 
a  floating  electrode  placed  between  each  transducer  electrode  (the 
split-isolated  transducer,  Fig.  12  inset).  The  relationships  between  the 
unconnected  electrode  voltages  (voltage  tap  weights)  are  not  known,  as  was 
the  case  for  the  transducer  discussed  in  section  III.,  therefore  the 
current  (ie.  admittance/voltage)  of  an  electrode  in  this  array  of  unknown 
voltages  cannot  be  calculated  directly  (as  was  the  case  in  (11)). 


TA 


.3 


j 


The  derivation  of  [M]  begins  by  cascading  single  electrode  matrices 
[m]  which  have  no  capacitive  electrostatic  admittance  (C=0).  This  general 
matrix  (similar  to  (22)) 
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models  all  the  relationships  (except  capacitive  coupling)  between  the 
electrodes.  The  admittance  terms  enclosed  by  the  dashed  line  are  the 
cascaded  complex  admittances  (G  +  JB  ) .  Added  to  this  sub-matrix  is  an 
inter-electrode  capacitive  admittance  (calculated  from  equations  (8),  (9) 
in  [2])  which  models  the  capacitive  coupling  between  the  electrodes. 


The  first  step  in  the  solution  of  the  floating  electrode  voltages  is 
to  combine  the  rows  and  columns  associated  with  the  transducer  (connected) 
electrodes.  This  is  done  using  the  same  technique  as  in  section  III.B. 

The  matrix  is  now  a  function  of  two  types  of  acoustically  induced 
potentials,  a  terminal  (transducer)  voltage  (VT)  and  a  set  of  floating 
electrode  voltages  (V?1  through  VPf ) . 
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(39) 


Note,  there  may  be  more  than  one  terminal  voltage  if  the  transducer  has 
multiple  groups  of  taps  connected  together  (ie.  group-type 
unidirectional) . 

An  examination  of  (39)  shows  that  the  induced  currents  (and  hence,  the 
floating  electrode  voltages)  are  functions  of  each  other,  the  terminal 
voltage  and  the  incident  acoustic  waves.  There  are  three  different  sets  of 
floating  voltages  corresponding  to  the  three  possible  excitation  sources 
(VT,  A1 ,  BN+1 ) ,  each  taken  one  at  a  time.  For  example,  the  value  of  is 
defined  as  the  traveling  acoustic  wave  emerging  from  the  left  port  due  to 
the  terminal  voltage  applied  with  the  incident  waves  absent. 

The  second  step  is  to  use  one  of  the  excitation  sources  as  the 
appropriate  bondary  condition  and  solve  the  admittance  sub-matrix  for  the 
three  sets  of  floating  voltages.  This  is  done  by  multiplying  the  first 
three  columns  (39)  by  the  appropriate  boundary  conditions  (each  0  or  1). 

The  Gauss-Seidel  iteration  method  [25]  is  used  three  separate  times  to 
solve  this  new  f  by  f+3  sub-matrix  for  the  unknown  floating  electrode 
voltages.  The  inversion  of  the  admittance  matrix  is  not  required  and 
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consequently  very  large  matrices  may  be  analyzed. 

The  third  step  is  to  solve  for  the  device  matrix  [M]  by  using  the 
floating  electrode  voltages  in  (40)  to  combine  the  columns  and  rows  as  in 
section  III.  (for  an  unconnected  electrode  b^  is  0).  This  is  done  for  each 
of  the  three  excitation  sources. 

As  an  example,  consider  the  solution  of  the  first  column  in  the  device 
matrix  (M^,  M21 ,  M^).  These  scattering  parameters  are  defined  as  the  two 
acoustic  waves  and  terminal  current  produced  by  an  incident  acoustic  wave 
from  the  left  when  the  device  is  shorted  (in  this  split-isolated  case,  a 
shorted  transducer)  and  the  incident  wave  from  the  right  is  not  present. 

Let  Bn+1  and  VT  equal  zero  and  the  incident  wave  (A1 )  equal  one. 
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After  multiplication  of  the  first  three  columns  (accounting  for  the 
boundary  conditions)  the  sub  matrix  is 
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After  solving  the  floating  voltages  from  these  f  equations,  the  desired 
coefficients  are  the  first  three  rows  (m=1,  2,  3)  multiplied  by  the 
external  boundary  conditions  and  the  potentials: 


MH  =  mn‘1  +  m12#0  +  mi3*°  +  mi4*VF1  + 
M21  =  m21#1  +  m22*°  +  m23*°  +  m24*VF1  + 
m31  =  *1  +  m32»0  +  “33*°  +  “34*^1  + 


+  m1f 

+  m2f 
+  m3f 
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(42) 


Similar  steps  are  taken  to  derive  the  other  six  device  matrix  [M] 
parameters.  This  matrix 
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completely  defines  the  transducer  with  floating  electrodes  and  is  used  to 
derive  the  electrically  loaded  parameters  via  (29-31). 


VI .  Conclusion 


This  work  has  modified  and  corrected  the  phased  periodic  array 
scattering  parameters  for  use  as  a  single  electrode  (unit  cell)  scattering 
matrix.  The  analytical  expressions  upon  which  this  matrix  is  based  are 
functions  of  material  constants,  electrode  geometry  and  frequency.  These 
single  electrode,  three  port  scattering  matrices  are  acoustically  cascaded 
to  produce  an  array  scattering  matrix  which  completely  describes  the 
device,  be  it  reflector,  transducer  or  combination.  The  analysis  takes 
into  account  the  effects  of  both  piezoelectric  and  mechanical  scattering 
and  provides  all  device  transfer  functions  versus  frequency.  By  virtue  of 
the  scattering  matrix  approach,  the  final  3-port  matrix  includes  all  finger 
interactions  such  as  acoustic  reflections  and  electrical  regeneration 
regardless  of  individual  electrode  circuitry.  The  calculation  is  a 
straight-forward  procedure  and  is  applicable  to  any  uniform  beam-width, 
periodic  electrode  configuration  (including  the  case  of  transducers  with 
floating  electrodes) .  Good  agreement  with  experimental  data  confirms  this 
approach  to  be  a  universal  and  accurate  analysis  tool. 


50 


Appendix  I. 


A.  Derivation  of  radiation  conductance  (G  ) 

3L 


The  conservation  of  energy  for  a  lossless  symmetrical  interdigital 
electrode  requires  the  electrical  input  power  to  equal  the  total  acoustic 


output  power.  That  is,  with  a  voltage  V  applied 

,2 


Ga  V*  =  Re(m33)  V*  _  1?^!^  #  !V  m^ 


,2 


2  Z, 


2  Z, 


(Al) 


dividing  both  sides  by  _vf  gives 

2 


Ga  =  m*i3  +  “23  m*23 


(A2) 


replacing  m*13  and  m*23  with 


°13  =  "Z0  m31  n23  =  "Z0  m32 


(A3) 


from  appendix  I.  B.  gives  equation  (7) 

Ga  s  “*5(n13m  31  +  “23“  32)' 


(AA) 


^  jk.ii  1  -  9- 
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B.  Derivation  of  the  relationship  between  and 

The  relationship  between  mixed  units  and  pure  units  scattering 
matrices  has  been  discussed  for  the  two  port  case  by  Viggh  [26].  This 
appendix  extends  his  work  to  the  three  port  case  in  order  to  derive  the 
relationship  between  the  mixed  units  m^  and  m^.  For  clarity  of 
presentation,  the  mixed  units  and  pure  acoustic  scattering  matrices  will  be 
designated  by  [m]  and  [s],  respectively. 


The  variables  of  the  mixed  units  electrical  port  are  voltage  and 
current.  Their  relationship  to  an  acoustic  representation  is  via 


V3  3  a3  +  b3 

z0  *3  =  a3  "  b3 


(A5) 


where  ZQ  ia  the  electro-acoustic  charactersitie  impedance  (8).  The  problem 
is  to  write  the  three  port  acoustic  scattering  matrix  (where  a^  and  bt  are 
the  input  and  output  acoustic  wave  potentials,  respectively) 

311  a12  313  a1 

S21  s22  S23  a2  I  (A6) 

S31  s32  *33 
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L  3  J 

in  terms  of  the  mixed  units  matrix  (similar  to  (1)) 
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Starting  with  the  third  nixed  units  equation,  insert  (A5)  for  the  current 


and  voltage 


a3  '  b3  8  “31a1  +  *32*2  +  n33(a3  *  V 


Rearranging  terms  and  solving  for  produces  the  third  acoustic  scattering 


equation 
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The  next  step  uses  the  first  mixed  equation  to  produce  the  first  acoustic 
equation 


b1  =  m11al  +  ml2a2  +  ®13v3 

8  m11a1  +  mi2a2  +  n13(a3  +  b3)- 


(A10) 


From  (A6)  we  know 


b3  8  a3lal  +  a32a2  +  s33a3* 


(All) 


and  therefore, 


b1  =  (m11  ♦  “13*31  )a1  +  (n12  +  m13S32)a2  +  m13°  +  s33)a3 


(A12) 


A  similar  equation  exists  for  bg. 
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In  summary,  the  acoustic  scattering  matrix  [s]  (A6)  in  terms  of  the 
mixed  units  representation  is 
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A  symmetrical,  lossless,  reciprocal  3-port  has  the  property  that  s^  = 
•  Therefore,  using  the  transformation  C  A 1 3 )  tie  find 
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producing  the  desired  result 


"l3 


(A15) 


it m  1 1 
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Appendix  II.  Scattering  Matrix  Cascading 

The  problem  addressed  in  this  appendix  is  to  find  the  resultant 
acoustic  potentials  B1f  a^  and  the  currents  produced  for  the  case  where  the 
incident  acoustic  potentials  f  Es  and  the  voltages  applied  are  known.  In 
this  case  the  imaginary  admittance  is  known  and  the  largest  matrix  is  N+2 


The  first  scattering  matrix  (which  could  be  the  entire  matrix  [M]  to 
the  left  of  this  point)  is  written  in  equation  form  as 


B1  =  M11A1  ♦  M12B2  *  M13VT 
A2  =  M21A1  +  M22B2  +  M23VT 
I1  8  M31A1  +  M32B2  +  M33VT 


(A16) 


where  the  have  been  multiplied  by  the  tap  weight  (ie.  the  voltage 
applied  as  a  fraction  of  the  terminal  voltage) .  The  set  of  equations 
characterizing  the  tap  to  be  appended  i3 

B2  8  m1 1 A2  +  m12B3  +  m13VT 

A^  =  ®2iA2  +  m22B3  +  m23VT  (A17) 

1Z  ~  nl3 1 A2  +  ®32B3  +  “33VT 


Comparing  the  matrix  equations  shows  that  there  are  two  equations  that 
contain  A.,  in  terms  of  Bg  and  the  known  quantities  (incident  acoustic  waves 
and  voltage  tap  weights).  These  are  solved  as  follows,  taking  the  second 


set  first 


=  ^  mi2B3  *  ®1  3^T 


inserting  an  equation  for  A^  gives 

®2  *  b11^M21A1  +  M22B2  +  M23^  +  b12B3  +  ®13VT 


(A  18a) 


(A 18b) 
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and 


B,,  -  mn^M21A1  +  M23^  +  m12B3  +  a13VT 

1  -  n11M22 

(A18c) 

which  is  written  as 

B2  s  Ciai  +  C2B3  +  C3VT 

(A19) 

where 

C,  =  m11M21 

1  "  m11M22 

(A20a) 

C„  3  n12 

(A20b) 

C_  =  ®11M2^  +  mnVT 

1  -  “1 1M22 

(A20c) 

Similarly 

*2  *  "21*,  *  "22B2  *  M23VI 

(A21a) 

is  written  as 

A2  =  DiAi  ♦  ^2®3  +  B3VT 

(A21b) 

where 

D1  3  M21  +  M22C1 

(A22a) 

°2  ■  "22c2 

(A22b) 

°3  *  M22C3  *  "23 

(A22c) 

These  solutions  are  inserted  into  the  proper  equations  in  the  two  matrices 
(A16,  A17 )  and  a  new  3  by  4  array  matrix  is  created  in  terms  of  the  new 
input  parameters 

B1  =M11A1  ♦M12(C1A1  ♦  C2B3  -  V  ♦  Mi3Vt 

Aj  s  m21(D1A1  +  B3D2  +  Dj)  ♦  ®22B3  +  m23VT  (A23) 

i,  =  m31a1  ♦  m32(c1a1  ♦  c2b3  ♦  c3)  +  m33vt 

I2  s  ®3^(D^A^  ♦  D2B3  +>  03)  +  1B32B3  +  ^33^T 
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Therefore  a  new  M  matrix  of  the  form 


M11  +  C1M12 

C2M12 

C3M12 

♦M13 

D1m21 

D2®21  +  m22 

D3m21 

+  m23 

M31  +  C1M32 

C2M32 

C3M32 

+  M33 

_fl“31 

D2m31  +  “32 

D3“31 

+  m33 

(A24) 


has  been  created  from  the  original  two  scattering  matrices.  This  matrix 
relates  the  input  acoustic  potentials  (A^B^)  and  the  terminal  voltage 
applied  (VT)  to  the  output  acoustic  potentials  (A^,  and  the  currents 
produced.  The  procedure  can  be  carried  out  for  N  taps  with  a  resulting  N+2 
by  3  scatter  matrix  for  the  entire  array. 
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